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FROM PARTICLE COUNTING TO GAUSSIAN TOMOGRAPHY 


K. R. PARTHASARATHY AND RITABRATA SENGUPTA 

Abstract. The momentum and position observables in an n-mode boson Fock space r(C"') 
have the whole real line K as their spectrum. But the total number operator N has a discrete 
spectrum Z_|_ = {0,1,2, - • • }. An n-mode Gaussian state in r(C") is completely determined 
by the mean values of momentum and position observables and their covariance matrix which 
together constitute a family of n(2n -I- 3) real parameters. Starting with N and its unitary 
conjugates by the Weyl displacement operators and operators from a representation of the 
symplectic group Sp{2n) in r(C"’) we construct n(2n -I- 3) observables with spectrum Z_|_ but 
whose expectation values in a Gaussian state determine all its mean and covariance parameters. 
Thus measurements of discrete-valued observables enable the tomography of the underlying 
Gaussian state and it can be done by using 5 one mode and 4 two mode Gaussian symplectic 
gates in single and pair mode wires of r(C"') = r(C)®”. Thus the tomography protocol admits 
a simple description in a language similar to circuits in quantum computation theory. Such 
a Gaussian tomography applied to outputs of a Gaussian channel with coherent input states 
permit a tomography of the channel parameters. However, in our procedure the number of 
counting measurements exceeds the number of channel parameters slightly. Presently, it is not 
clear whether a more efficient method exists for reducing this tomographic complexity. 

As a byproduct of our approach an elementary derivation of the probability generating 
function of in a Gaussian state is given. In many cases the distribution turns out to be 
infinitely divisible and its underlying Levy measure can be obtained. However, we are unable 
to derive the exact distribution in all cases. Whether this property of infinite divisibility holds 
in general is left as an open problem. 


1. Introduction 

It is in the nature of quantum theory that properties of the state of quantum systems can be 
inferred from measurements of observables of physical signihcance taking values in a discrete set 
or a continuum. From an experimental point of view it is natural to seek as much information 
as possible from the discrete measurements. A typical measurement of the discrete type is 
counting the number of particles of a particular type. Suppose the unknown state of a system 
can be described in terms of some parameters which constitute a manifold of dimension k. 
Then it is natural to look for k discrete-valued observables from whose expectation values one 
can determine the values of these parameters. Of course, such observables should be physically 
meaningful and also experimentally measurable. This has been extensively studied in the book 
[PROdj . 

In this article we explore this problem of determining the state when it is known that the state 
is Gaussian. Consider the Hilbert space or equivalently, the boson Fock space F(C"') 

over the n-dimensional complex Hilbert space C”. A Gaussian state of n-modes in is 

completely described by its momentum and position mean values and a 2n x 2n covariance 
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matrix. Thus, an n-mode Gaussian state is determined by n{2n + 3) parameters. Here one has 
observables a^aj, the number of particles in the j-th mode for each j = 1,2, ■ ■ ■ ,n and also their 
unitary equivalents in different frames which are obtained by Weyl (displacement) operators as 
well as unitary operators which implement the symplectic linear transformations in the position 
and momentum observables obeying the canonical commutation relations. Using these resources 
we shall construct n{2n + 3) number observables which have the discrete spectrum {0,1, 2, • • • } 
and the property that all the means and covariances of the unknown Gaussian state can be 
easily determined from their expectation values. Measurements on these number observables on 
an ensemble of such a Gaussian state and the law of large numbers can be used to estimate the 
unknown parameters. In the process of such an investigation we shall determine the probability 
generating function of the distribution of the total number observable its mean and 

variance in any hxed Gaussian state. 

Following Heinosaari et al |HHWin] a Gaussian channel with n degrees of freedom is deter¬ 
mined by a pair {A, B) where H is a 2n x 2?7, real matrix, B is a. 2nx2n real positive semidehnite 
matrix satisfying the matrix inequality 


B -|- J2nA — J2n) ^ 0, 


where J 2 n is dehned in equation (2.12). Thus such a channel is determined by 4n^-|-n(2?7,-|-l) real 
parameters. Such a channel yields an output Gaussian state for any given input Gaussian state. 
By choosing a few appropriate input coherent states and performing a Gaussian tomography 
on the output Gaussian states, we show how the matrices A and B of the quasifree channel can 
be estimated. We use coherent states as they constitute an important class of mathematical 
objects, which are easy to realize experimentally. Further, creation of different coherent states 
with different amplitudes and phases can be possible from the same experimental set-up. We 
note that, a similar approach of using coherent states for tomography (although in a different 
scheme) has also been taken by Lobino et al |LKK'*~n8] , and was further developed in jR.KSM 
IWYH+I3j . 


T1 


Bosonic particle counting is an important tool in quantum optics, both for theory and for 
experiments |SW87bllSW87aj . Gaussian states, channels and their applications have been used 
extensively in quantum information theory. These concepts have been studied in detail in the 
book of Holevo |Holl2j and in the survey article by Weedbrook et al |WPGP'*~12j . We also 
refer to the books by de Gosson [dGOb] and Parthasarathy |Par92j for the connection between 
symplectic geometry and quantum stochastic calculus. We refer to the survey article by Lvovsky 
and Raymer |LR,n9j and the references therein, for experimental processes of continuous variable 
tomography. 

We organise the paper as follows. To increase the readability of the article, in we write 
a short introduction to notions like exponential vector, Weyl operator, Gaussian state, Fourier 
transform, Gaussian channels and other necessary concepts which will be used in subsequent 
sections. In this venture, we mostly follow the approach taken in the papers |ADMS9^ IParlO( 
IParl4aj . In (j^we derive the basic formulae for expectation and variance of the total number 
operator which will be used in (Q to estimate an unknown Gaussian state. In we use the 
tomographic method derived for Gaussian states in (Qto estimate the unknown parameters of 
a Gaussian channel. 
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2. Notation and preliminaries 


In this section we give a short survey on Gaussian state and other necessary concepts. Though 
these have been extensively studied in various references, we follow the method and notation 
adopted in the book |Par92] and the papers |Parini IParl3| IParl4bl IParldaj . 


2.1. Exponential vector. Let H be a finite dimensional complex Hilbert space. When dimH = 
n and H is identihed with C” we express its elements as column vectors z = ^2, ■ ■ • ! ZnY 

with Zj being complex scalars and the scalar product between two elements z and z' as 

(z|z') = 

i=i 

Dehne the boson Fock space T(FL) over FL by 

(2.1) r(H) = c©H©H®^©---©H®"©--- 

where (s)^ denotes r-fold symmetric tensor product. Elements of the subspace 'H®'" in T(FL) are 
called r-particle vectors and elements of the form 

Mo © © • • • © Mr © • • • 

where all but a finite number of Mr’s are null, are called finite particle vectors. Finite parti¬ 
cle vectors constitute a dense linear manifold T in T{'H). For any u E FL we associate the 


exponential vector e(u) 

in r('H) defined by 


(2.2) 

U®2 

e(u) = 1 © u © —= © • • • 

© ^ © • • 

Vr! 

Then 



(2.3) 

(e(u) e(v)) = exp(u v) 

Vu, w E Ft 


The linear manifold S generated by all the exponential vectors is called exponential domain. 
The two dense linear manifolds £ and iF are useful domains for constructing several operators 
of physical significance. 


2.2. Weyl operator. For any u G "H we associate the Weyl displacement operator hP(u) by 
putting 

(2.4) iy(u)e(v) = + v) 

for all \ E FL, observing that iy(u) is scalar product preserving on £ and therefore extends 
naturally to r('H). The Weyl operators obey the multiplication property 

(2.5) W(u)W(v) = + v) 

for all u, \ eFL and yield a strongly continuous, irreducible, factorizable and projective unitary 
representation of the additive group Ft. By ‘factorizable’ we mean the property that under the 
isomorphism between T{FLi ®FL2) and r('Hi) ®V{FL2) through the identification Xe(ui ©U2) = 
e(ui) © e(u2) one has 


XW(U1 © U2)X-^ = W(ui) © W(U2). 
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If u —)■ VT (u) is another str ong ly continuous map from Tf into the unitary group of a Hilbert 
space /C such that equation (2.5) holds with W replaced by W and W{-) is irreducible, then 
there exists a unitary isomorphism V : r('H) —)■ /C such that 


VW{u)V-^ = IH(u) Vu G H. 

Thus the Weyl operators constitute a unique multiplicative family up to unitary equivalence 
but with the presence of the factor exp — ?Im(u|v) in equation (2.5). We call u —)■ iy(u) the 
Weyl representation of H and we shall exploit its properties to define a natural quantum Fourier 
transform for states in T{'H). For now, we shall introduce some basic observations arising from 
the Weyl operators. 

For any hxed u G Tf, the map t W (fu) yields a strongly continuous one parameter unitary 
group as t varies in M. By Stone’s theorem |Par92] there exists a self adjoint operator p(u) such 
that 


(2.6) W{tu) = t G M. 

Dehne the operators 

g(u) = -p{tu), 

«(u) = -(g(u)+?p(u)), 

= ^(g(u)-zp(u)). 

All these operators have domains including the exponential domain £ and the domain if of 
finite particle vectors. Indeed, any finite linear combination of these operators have the same 
property and we denote their respective closures by the same symbols. With this convention 
one has 

W(u) = 

When "H = C” = M” + and u = x + ?y with x = Reu, y = Imu we also have 

W(u) = W(x + ty) = 

Furthermore, one has the following commutation relations. 


[P(u),p(v)] 

[?(u),g(v)] 

[g(u),p(v)] 

[a(u),a(v)] 

[a^(u),a^(v)] 

[a(u),a^(v)] 


2 z Im(u|v), 
2 t Im(u|v), 
2i Re(u|v), 
0 , 

0 , 

(u|v) 


on the domains £ and if for all u, v in "H. 
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Choose and fix an orthonormal basis {ej}, j = 1, 2, • ■ ■ in "H, and define 
Pi = 

% = ^ ^ ^ ^ ■ 

Then one has the canonical commutation relations in the form 

\Pr,Ps] = [gr, Qs] = 0 , [qr,Ps\ = l^rs, 

or equivalently, 

[(X,’, II5] O5 [Hr? II5] 

in the domains S and J^. The observables Pi,P 2 ,-'' cire called momentum operators and 
qi, g 2 , • ■ ■ are called position operators in the basis {ej, j = 1, 2, ■ ■ ■ }. 


2.3. Quantum Fourier transform. For any trace class operator p in r('H) its quantum 
Fourier transform or simply Fourier transform p on "H is defined by 

p(u) = TrphF(u), ueFL. 

Then p is a bounded continuous function of u satisfying p(0) = Trp. If p is positive then p 
obeys the Bochner property: for any finite set {cr, r = 1,2, ,k} of scalars and elements 

{ur, r = 1, 2, • • ■ ,k} in FL one has the inequality 

CrCs exp(xIm(Ur|Us))p(Us — Ur) > 0. 
r,s 

Conversely, if (p is a continuous function with <p(0) = 1 and <p satisfies the Bochner property 
above then there exists a unique state p such that p = ip. When Ft = one has the Fourier 
inversion formula: 

^ J du 

where du denotes integration with respect to the 2n dimensional Lebesgue measure in with 
du = dxdy, u = x + xy, x = Re(u), y = Im(u). 

With the help of Fourier transform we shall now construct a natural Hilbert space isomor¬ 
phism between the Hilbert space of Hilbert-Schmidt operators in F('H) and the Hilbert space 

L2(M2n^_ 

Proposition 2 . 1 . Let FL = C”. Then 

— j exp [—||w|p -I- (u|w) + (w|v)] dw = exp(u|v). 

Proof. Immediate from standard formulae for Gaussian integrals. □ 

Proposition 2 . 2 . Denote by F^iFL), the Hilbert space of square integrable functions on Ft with 
the scalar product 

(/|g) = j 
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and by 82 ( 1 '(Ti)) the Hilbert space of all Hilbert-Schmidt operators on T(H) with the scalar 
product 

(P1IP2) = Trplps- 

Then there exists a unique Hilbert space isomorphism F : B 2 (J'{'H)) —)■ L‘^{fH) such that for any 

u, \ eH, 

(2.7) (F(|e(u))(e(v)|)) (w) = Tr|e(u))(e(v)|PF(w) Vw G 77. 


Proof. The right hand side of equation ( |2.7 ) is equal to 

(e(v)|I.r(w)|e(u)) 


(2.8) 


= e 


= exp 


e(v)|e-<'^l’">|e(u + w)) 

w|u) + (v|u) 


■^||w|p + (v|w) 


If we put pj = |e(uj))(e(vj)|, j = 1,2 then 
(2.9) 


Ttp[p 2 = exp [(ui|u 2 ) + (v 2 |vi)]. 

On the other hand the scalar product between Tr|e(uj))(e(vj)|hh(w), j = 1,2 reduces by 
Proposition 2.1 and equation (2.8) to the right hand side of (2.9). Now we observe that rank 
one operators of the form |e(u))(e(v)|, u, v G 77 constitute a total set in i32(r(77)). Thus F 
dehned by (2.7) extends uniquely to t he w hole of i32(r(77)). On the other hand functions of w 
of the form on the right hand side of (2.8) constitute a total set in L^(77). Thus F extends to 
a Hilbert space isomorphism between H2(r(77)) and □ 

2.4. Gaussian state. 

Definition 2.1. A state p in r(77) with PL = C"' is called an n-mode Gaussian state if its 
Fourier transform p is given by 

T 


( 2 . 10 ) 


p(x + zy) = exp 


—z-\/2(l^x — m'^y) — 




for all X, y G M” where 1, m are elements of M"' and S is a real 2n x 2n symmetric matrix 
satisfying the matrix inequality 

( 2 . 11 ) 2 ^ + zJ 2„>0 

with 


( 2 . 12 ) 


J2n, 


0 

In 


-In 

0 


In being the identity matrix of order n. 


Remark 2.1. Equations (2.10)-(2.12) have been written keeping in mind the orders of the 
canonical momentum and position observables as pi, p 2 , ■ ■ ■ , Pn, <? 2 , ■ ■ ■ , Qn- Sometimes it 
is more convenient to distinguish the different modes of a Gaussian state by using the order 
P 15 P 2 , Q' 2 , ■ ■ ■ ^ Pm An- This is usually achieved by employing the permutation 

_ A 2 3 4 ■ ■ ■ 2n - 1 2n 

\1 n 2 n +1 ■■■ n n + n, 
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Then the right hand sides of (2.10 )-(2.12) are obtained by changing (x^, y^), (1^, m^) and S re¬ 
spectively to {xi,yi,X2,y2, - ■ ■ ,Xn,yn), (^1, , ^2, " 12 ,''' , 4 , "In), and with J2n replaced 

by 

0 -1 
1 0 

0 -1 

J 2 n = 1 0 

0 -1 
1 0 

By abuse of notation we may denote both J2n and J2n by the same symbol J2n- 

We choose the canonical orthonormal basis ej = (0, 0, • • ■ , 1, 0, ■ ■ ■ , 0)^ with 1 in the j-th 
position for j = 1, 2, • • •, then the momentum and position operators satisfy the relations 

Trpjp = Ij, Tr qjp = rrij 

and S is the covariance matrix of {pi,P2, ■ ■ ■ ,Pn, —Qi, —<12, ■ ■ ■ , —Qn) in the state p satisfying 


(2.10 2.12). Whenever 2.10 is satisfied we write 

p = Pg(l,m; ^). 

Thus p is completely described by n{2n + 3) parameters. 

Proposition 2.3. Let pg(lj, mj; S'j), j = 1,2 he two n-mode Gaussian states. Then 

exp 


-1 

1 

1 

T 

(S + T)-^ 

1- 

1 - 

1 

1 _ 

Zr 

-l^mi - m2j 


-l^mi - m2jj 


(2.13) Tr pg(li,mi;^)pg(l 2 ,m 2 ;T) = - 

A/det(S T) 

Proof. Any state in r(C”) is a positive operator of unit trace and hence a Hilbert-Schmidt 


operator. Thus by Proposition 2^ we have 
Tr pg(li, mi; S)pg{\ 2 , m 2 ; T) 

1 


exp 


7\/2((1i - l 2 )^x - (mi - m 2 )'^y) 


{S + T) 


dx dy. 


The rest follows from the standard formula for the characteristic function of a multivariate 
normal density function in statistics. □ 

Proposition 2.4. For any u G C” with x = Re(u), y = Im(n) 

IV(u)pg(l, m; S') W(u)t = pg(l', m'; S) 

where 

h = 1 + V^y, m' = m -|- a/2x. 

Proof. Immediate from Corollary 3.3 in |Parinj . □ 























K. R. PARTHASARATHY AND RITABRATA SENGUPTA 


We denote by Sp{2n) the symplectic group Sp{2n, M) of all real 2n x 2n matrices L satisfying 
the relation 

L J2nL •J2n- 

Let r(L) be the unitary operator in r(C’^) which is unique upto a scalar multiple of modulus 
unity and satishes the relation 

r(L)iy(x + zy)r(L)-^ = W{x.' + ?y'), Vx, y G M”, 

where 

' ' 'x' 

y' 

If U is any n x n unitary matrix and U = A + iB where A = Re(17) and B = lm{U) then 
the 2 n x 2 n matrix 

'A -Bl 
B A 


L = 


is an orthogonal matrix which is also an element of Sp{2n). We denote the corresponding r(L) 
by r(f/) and call it the second quantization of U. We can realize r(f/) as the unique unitary 
operator satisfying 

r(t/)e(u) = e{Uu) Vu e C”. 

With these notations we have 
Proposition 2.5. For any L G Sp{2n) 

r(L)p,(l,miS)r(L)t=p,(l',m'iS') 

where 


1 ' 


-m 


r 


-m 


= {L-r 

s' = {L-^fSL-\ 

Proof. This is Corollary 3.5 of [Par M- 

3. Particle counts and their statistics in a Gaussian state 


□ 


Let pg(l, m; S) be an n-mode Gaussian state in r(C"') and whose Fourier transform is given 
by equation (2.10). Dehne the observables 

Nj = a]aj = ]^{p] + - 1), 1 < j < n. 

n 

N = Ey- 

i=i 

Both Nj and N are observables with spectrum {0,1, 2, ■ ■ ■ }. Nj is called the number operator 
or observable which counts the number of particles (photons) in the j-th mode. Since the 
Nj^s commute with each other they have a joint distribution in the state Pg(l, m;S') with 
support in {0,1, 2, • • ■ Using Proposition 2.3 we shall derive a formula for the probability 
generating function of this joint distribution and arrive at some natural corollaries. To this end 
we introduce the Gaussian state 
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where 


p(0,0;T) = JJ(1 - 
i=i 


T 


"1 f l+e-*l A 

2 J 


h 



/1 + e h \ 
V1 - e“h j 


1 

2 


/ l+e~*” \ 
1—e“*n j 


h 


h, l<j<n 


tj > 0 Vj, 


with D indicating the diagonal block matrix with blocks of order 2x2 and the diagonal entries 
being enumerated within ( ). These are the well-known thermal states. It follows immediately 

from Proposition 2.3| equation(|2.13[) that 


exp 


Trpg(l,m; S)e 




-m 


S + D 


l<J<n 


-1 


-m 


m=i(i 


e 


^ + l<j<n 


det 


Substituting Nj = Xj = e b we get 

(3.1) 

exp 

Trpg(l, m; S)x^^ ■ ■ ■ = - 



S + D 


l-Xj 


det 


S + Dil 


h, l<j<n)) 


for 0 < Xj < 1 Vj. The right hand side of this equation is nothing but the probability generating 
function of the joint distribution of the observables Aj-, 1 < j < n. 


Theorem 3.1. Let pg(l,m; S) be an n-mode Gaussian state in r(C”') whose covariance matrix 
S has eigenvalues Ai > A 2 > ■ ■ ■ > A 2 n with respective real eigenvectors bi,b 2 ,-- - , b 2 n and 


momentum position mean 


1 

m 


satisfying 


Suppose 



2 n 




aj = 



Xj -f- 


1 ’ 

2 


1 < j < n. 
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Then the probability generating function Gn{x) of the distribution of the total number operator 
N in the state Pg(\, m; S) is given by 

Gn{x) = Tt pg{\,ni-,S)x^ 

2n 


(3.2) 



1 -g; 

1 — ajX 


exp 


1 ^ 2 ( 1 -«i) 
2 ^ 


(1 - ajx) 


0 < a; < 1. 


Proof. Putting Xj = x for all j in (3.1) and making use of the eigenbasis and eigenvalues of S 
we see that (3.1) becomes 

, -1' 


2n 


exp 




Gn{x) - P[ . - . 

i=i Y (1 ~ + 2 ^) 

The rest is elementary algebra using the dehnitions of aj in terms of \j for every j. 


□ 


Corollary 3.1. 
(3.3) (» 


The probability distribution of N in the state pg{\, m; S) satisfies the following: 

2n 


(3.4) 


(3.5) 


Ui 


[III] 


Pt{N = 0) = P[(l 
1=1 




aj) exp 




a,- 


/ 1\ 


2 " 

Tr 5 - - + 


[ V 2; 

{-mj 



Variance{N) = -^Tr (“S' — ^ 


-2 


-m 




-m 


Proo f. Property (i) follows by putting a: = 0 in (3.2). Properties (ii) and (iii) are obtained from 
(3.2) by taking the logarithm of G^, differentiating twice and taking lim 3 ._^i ^(logG(a;)) and 

□ 


lim^^i ^(logG(a:)). 


Remark 3.1. Equations (3.3)-(3.5) show that the probability of presence of a particle and 


the expectation and variance of the total number of particles get enhanced when the Gaussian 
state has a nonzero momentum position mean vector. Indeed, the mean and variance of N tend 
to inhnity as the length of the momentum position mean vector increases to inhnity. Equations 


(3.3) and (3.5) indicate the possibility of estimating the mean and covariance parameters of 


a Gaussian state by measuring the number operator under different displacements. We shall 
discuss this approach to the tomography of a Gaussian state in great detail in the next section 


It may be noted that the parameters a,- in Theorem 3.1 satisfy the inequality |ad < 1 for 


every j. If 1 and m are null-vectors and aj > 0 for every j then Gn{x) assumes the form 

2n , 

G„(u=n()^“' 


1=1 


ajx 


and the corresponding distribution of iV is a convolution of 2n negative binomial distributions 
of index some of which may be degenerate at 0. In particular, it is inhnitely divisible. In 
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this case 

Variance(A^) — (N) = ^Tr ^‘5' — >0 

and the distribution exhibits a super Poissonian property. When S = \l 2 n the state becomes 
the vacuum state and N has the degenerate distribution, degenerate at 0. 

Now we shall analyse the distribution of At in a pure Gaussian state. In this case S = \L?"L 
for some element L of the group Sp{2n) and its eigenvalues can be expressed as 

2 ’ 2c 

where Cj > 1 for 1 < j < k. Then 

Oi2j-l = 
a2j = 

Or = 


i’ 2’ 

2C2’ 

’ 2 ’ 

2cfc’2’ 

’ 2 J 

Cj-l 

Cj+1 

> 0 

for 

3<j<k, 


1 —Cj 

1+Cj 

< 0 

for 

3<J<k, 


0 


for 

2k + 3 < r 

< 2n. 



Write 13j = a 2 j-i, 3 < j < k so that a 2 j = —(3j, 3 < j < k. Then the probability generating 
function of the total number operator N in (3.2) assumes the form 


(3.7) 


G„(i) = 


where 


(3.8) 

Gi{x) 

k 

/ 1-/^1 

1 - /3^x^ 



(3.9) 

G2 {x) 

1 ^ 

i=i 

7 

1_1 

(x- 1)(1-/3j) ^^{x-3){3 + I3j) 

3 — 13 jX 1 + 13 jX 




'2n 



(3.10) 

Gsix) 

= exp 

- y 

2 ^ ^ 

(x- 1) 



. j=2k+l 


where 0 < Pj < 3 for 3 < j < k. 
Writing 


(3.11) 

(3.12) 


= T|,.-i(l-ft), 

= ^2,(1+ft), 1 <!<*:, 


one can express G 2 {x) in (3.9) as 


(3.13) 


Goix) = 


1 - 1) + [7i(l - Pj) + + f3j)] {x - 1) 


^ Pj\lj 
exp - ^ -L^ 


j=i 


1 — /3'^x‘^ 


where 0 < /3j < 1, 7 j > 0, 6 j > 0 for 1 < j < fc. 

From (3.7) it follows that Gi(x) is the probability generating function of a convolution of 
probability distributions pj, 3 < j < k where the probability generating function of pj is equal 
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to 


(1 - 13 ]) 


oo 

r=0 


1.3.5.---(2r + l) 


T ! 




2r 2r 

j ’ 


l<j<k. 


In particular fij is an infinitely divisible distribution with support in {0,2,4, ■■■}. Equation 
(3.10) shows that 6 * 3 ( 0 ;) is the probability generating function of a Poisson distribution with 
mean value ^ Y],l =2 


=2A:+1 6- ■ 

If 7 j > Sj for every j = 1 , 2 , • ■ ■ ,k then ( 62 ( 0 ;) is clearly the probability generating function of 
an inhnitely divisible distribution with support in {0,1, 2, • • • }. Its Levy measure can be easily 
read off from (3.13). Under this assumption, i.e. jj > 6j for each 1 < j < k it follows that N 
has an inhnitely divisible distribution in the pure Gaussian state we started with. 

However, we do not know the answer to the question whether the distribution of the total 
number operator N in every Gaussian state is inhnitely divisible and hence of a mixed Poisson 
type. 


4. From particle counting to the tomography oe a Gaussian state 


A Gaussian state with n modes can be constructed if its momentum and position means 1, 
m and its covariance matrix S are known. Our aim is to express these parameters in terms of 
the expectation values of conjugates of the total number operator by a few elementary gates 
in the Hilbert space P(C’^). To this end we shall make use of the Weyl displacement operators 
lU(u), M G C"" and the Gaussian symmetries r(L), L G Sp{2n) described in Section We 
start with an elementary but basic result for achieving this tomography of a Gaussian state. 
For any observable X denote by {X) its mean value in the state relevant to the context. 


Theorem 4.1. Let Pg{l,ni]S) be a Gaussian state and let N be the number operator in the 
n-mode Hilbert space r(C"). Then the following hold: 

(i) For all x, y G M"" 


(4.1) 


(lU(x + + W)> - (N) = l|xf + ||yf + V2{y^\ + x^m) 


(a) For any L G Spifln) 

(4.2) 

(r(L)'Afr(L)) -{N) = \ 




+ 


-m 


- hn 


-m 


Proof. We have from Proposition 2.4 and equation (3.4) 


(lU(x + *y)'*'A^lU(x + *y)) = Tr pg(l, m; S')lU(x + *y)'^iVlU(x + *y) 

= Tr lU(x + *y)pg(l, m; S')lU(x + *y)'^iV 


1 

2 


Tr f S' - ^ j + ||1 + V^yll^ + l|rn + V^x| 


Now subtracting the value of (N) given by (3.4) we obtain equation (4.1). 
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Similarly , we have from Proposition 2.5 and equation (3.4) 


Trpg(l,m;^)r(L)tiVr(L) = Ti T{L)pgil,ui; S)T{L)^N 

1 
2 


Tr 


L-^^SL-^ _ i ) + 


1 

—m 




1 

—m 


Now, subtracting the value of {N) given by (3.4) get get (4.2). 


□ 


For any unitary operator U in the 1-mode Hilbert space r(C) we say that the operator 

/ with U in j-th position, acting in the n-mode Hilbert space 
r(C”) = r(C) (g) r(C) 0 • • • (g r(C) is the gate U applied on the j-th mode and denote it by 

'-V-^ 

n times 

In Figure 0 we represent by following the notion of circuit diagrams in quantum 
computation. 


[/(i) = 



Figure 1. 

Here each wire stands for a F(C) and in the j-th wire the unitary operator U is applied. 

Similarly, if H is a unitary operator in the 2-mode Hilbert space F(C^) = F(C) 0 F(C) we 
construct the unitary operator in the 2 modes representing the Fth and the j-th wire. 

For example is represented in Figure 


H(1.2) = 



Figure 2. 

When i < j are not successive we can apply a permutation to make them successive, apply V 
and follow by the reverse permutation. One may use Figure]^ For achieving the tomography, 
we shall use only one and two mode gates. 

To begin with we consider the two 1-mode Weyl displacement operators W ( 2 “ 2 *) and W ( 2 “ 2 ) 
where t = \/—1. Put 

(4.3) Gp = hF(2-0) 

(4.4) G, = hF(2-i) 

If = (0, • • • , 0,1, 0 ■ ■ ■ , 0)^ with 1 in the j-th position then applying Gp and Gq in the j-th 
mode is equivalent to using the displacement operator W{2~Hej) and W{2~^ej) respectively. 
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i 

J 



Figure 3. 


Then equation (4.1) in Theorem 4.1 reduces to 


(4.5) I, = (gW'atgU)) - (AT) - 1, 

(4.6) rrij = - {N) - 1. 


Furthermore (3.4) implies 

(4.7) 


Tr S = 2(Af) 


m + n. 


In other words, the measurement of counting observables iV, Gp'^ NGp \ G^/^ NG^q\ 1 < j < n 
which constitute a set of cardinality 2n + 1 yields the 2n + 1 parameters Ij, rrij (1 < j < n) 
and Tr S concerning the Gaussian state pg(l,m; S). 

Till now, the only resources we have used are identity (or zero mode) gate and the one mode 
gates Gp and Gg. Now we shall pass on to 1-mode Gaussian symmetry gates. 

We consider the element L{x,a) G Sp{2) dehned by 


(4.8) 


L{x, a) 


cos a — sin a 


X 0 


cos a sin a 

sin a cos a 


0 i 

L xA 


— sin a cos a 


and the unitary operator 


X > 1, 0 < a < 271 


(4.9) Gsp{x,a) = r(r(L(a;,Q!))) 

where r(L) = (L“^)^ in any symplectic group. We now view Gsp{x,a) as a 1-mode gate and 
apply it in different modes. We express it by Figure where the box is in the j-th wire. 












Figure 4. 

Expressing the covariance matrix S of the n-mode Gaussian state as a block matrix 

bJ e {1,2, ••• ,n} 
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where each Sij is a 2 x 2 matrix we observe that Sjj is the covariance matrix of the j-th marginal 
Ganssian state and 


(4.10) 




% 


qt c. . 


Z < ] 


is the covariance matrix of the i j-marginal Gaussian state of Pg(l, m; S). 
Let 


(4.11) 


Sii = 


PP 

qp 


CP, 


pq 


U, 




with aqp = ^q. Thus Sjj has three parameters. Applying the gate Gip{x,a) dehned by (4.9) 
and Figure |4, and using part (ii) of Theorem 4.1 we get 

(4.12) 


{x^ cos"^ a + X ^ sin^ a — l)app + (x^ sim a + x ^ cos"^ a — l)<Jqq + 2(x 


X ^) sin a cos a 


= 2 [(^G^JJ{x, aYNG^JJ{x, a)) — (N)) + [1 — (x^ cos^ a + x sin^ a)]lj 


+ [1 — {x^ sin"^ a + X ^ cos"^ a)]m^ + 2ljmj{x 


X ^) sin a cos a 


Ghoosing a = 0, (4.12) becomes 


(4.13) 


{x^ - l)app + (x ^ - l)aqq 

= 2 ((G(^p)(x,0)^A^G(^p)(x,0)) - (N)) + (1 - x‘^)l‘^j + (1 - x-^)m 


Ghoosing x = a/2 and x = a/3 successively (4.13) yields two linearly independent equations for 
determining cXpp and aqq in terms of /•, rrij, {N) and (Gip^(x, 0)lA^Gip (x, 0) 


Now we go back to the equation (4.12) and choose x = a/2, a = j. Then we get 


^ WPP 


(J, 


qqJ 


-CP, 


pq 


= 2 


Ff i) > - w 


TT 


(4.14) 




Since (T„p, aqq and L, rrij have already been determined, (4.14) determines apq by using values 


of (TV) and (^G%^ (x/2, f)^iVGii^ {V2, f)). Thus Sjj can be completely determined by mea¬ 
surements of N using the gates Gsp{\/2, 0), Gsp(\/3, 0) and Gsp(V2, |) in different modes after 
knowing the vectors 1 and m. However, after determining Sjj for 1 < j < n — 1, in order to 
determine Snn h is enough to use only the two gates Gsp{y/2, 0) and Gsp{\/2, j) in the n-th 
mode because we already know Tr S form (4.7). Thus we need only (3n — 1) new measurements 
to determine the 3n parameters occurring in the block diagonals Sjj, 1 < j < n. 

Now it remains to determine for any i < j the off-diagonal block Stj. To achieve this goal we 
shall use a 2-mode gate of the form r(L), L G 5*^(4). We start with a unitary matrix of order 
2 of the form 


U = 


( 


1-/T a) ’ 

\a\ 
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where a = ai + ia 2 i /3 = (3i + i(32 with f3j being real. If we view f/ as a real linear 
transformation of we get an element of 5^(4) of the form 

oti —(y.2 /?! ~/^2 

a2 Oil (32 I3i 
—/3i —/?2 Ctl Cl2 

(32 ~(3l —Ot2 Oti 

where (9 is a real orthogonal matrix. Define 


(4.15) 


O = 


(4.16) 

We write 

(4.17) 


L{U, Xi,X2) = O 


Xi 




X2 


X2^ 


O^, Xi >1, X2 > 1. 


L{U,Xi,X2)^L{U,Xi,X2) = 


A B^' 
B C 


and note that A, C are 2x2 positive definite matrices and B is given by 

-(3iai{xi - X 2 ) + (320i2{x'l^ - X2^) -I3ia2{xi - x^^) - (32ai{x((^ - X 2 ) 
(32ai{xi - X2^) + (3ia2{x((^ - X 2 ) (320t2{xi - X 2 ) - (3iai{x((^ - x^^) 


(4.18) 


B = 


Using (4.10) and (4.17) we observe that 


Tr 


Sii Sii 


Sji Sjj 




'A B^' 



B G 


L{U,Xi,X 2 )'^L{U,Xi,X 2 ) = Tr 

(4.19) = Tr(5ii7l + SjjC) + 2Tr SijB 

where the first sum on the right hand side depends only on Su and Sjj which have already been 
determined in terms of the expectations of N and its conjugates by chosen one mode gates. In 
order to determine Stj we use the 2-mode gate 

(4.20) Gsp { U , xi , X 2 ) = r(r(L(U,xi,X2))) 


in the (f, j)-modes and use part (ii) of Theorem 4.1 to obtain the relation 


{G^’^{U,x„X2yNG^^^{U,x„X2)) - {N) 


1 

2 


Tr 


'Su 


-1 





A-I 2 B^ 

B G-h 


f h \ 

—mi 

k 

\-m 2 J 


71 - I 2 B^ 

B G-h 



/ \ 



—mi 



h 



\-rrL 2 j 



Using (4.19) this reduces to 

TrSyB = 

-tlV |S„(/1 -/j) + S^^(C - «| 

= f(U,Xi,X 2 ) say. 


(4.21) 
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When li, rrii, Ij, mj, Sa and Sjj are already determined, the term f{U, Xi,X 2 ) depends only 
on U, Xi, X 2 . Let 




7ii 7 i2 
721 722 


We now make fonr special cases for {U, Xi, X 2 ) 

1 1 
-1 1 


= H = 


(4.22) 


, Xi = 1, X 2 = 2. Pnt ri = f{H, 1,2). Then (|4.2l|) becomes 
1 1 


27ii “ ^722 — Li. 


(ii) U = H, xi = 1, X 2 = 3. Pnt r 2 = /(i/,1,3). Then (4.21) becomes 
(4.23) ' 

(hi) [/ = A- = i 


7ll ~ ^722 — ^ 2 - 


I 1 

— 1 —I 


, xi = 1, X 2 = 2. Pnt ra = /(iP, 1,2). Then (4.21) becomes 


(4.24) 


-(721 + 2712) = rs. 


(iv) U = K, xi = 1, X 2 = 3. Pnt = f{K, 1,3). Then (4.21) becomes 
(4.25) - ^(721 + 3712 ) = r4. 


The fonr eqnations (4.22)-(4.25) in the nnknowns 711 , 722 , 712 , 721 are linear and linearly 
independent. Thns they determine the matrix Sij for any hxed i, j. For this pnrpose we have 
nsed exactly fonr measnrements described by the fonr 2-mode gates for fonr parameters. 

In all we have nsed exactly (2n -|- 1) -|- (3 ?t. — 1) -|- = n(2n + 3) measnrements to 

determine the n{2n + 3) parameters of the Ganssian state Pg(l, m; S). 


5. Tomography of Gaussian channels 


An n-mode Ganssian channel /C(A, B) is described by a pair of real 2n x 2n matrices (A, B) 
where B is positive semidefinite and the following matrix ineqnality holds: 


B -|- J2nA — J2n) ^ 0 


with J 2 n as in eqnation (2.12). Thns }C{A,B) is determined by 6n^ -|- n real parameters. Snch 
a channel has the property that for any Ganssian inpnt state ^^(l,!!!;^), the corresponding 
ontpnt is again Ganssian and has the form Pg(}', m'; S') where 


(5.1) 

(5.2) 



A' = A'^SA+^B. 
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our aim is to determine A and B by performing tomography on the output states for a small 
number of coherent input statea |'0(u)) where 


(u))(^(u)| = pg y2x; 


where x = Re(u), y = Im(u). By (5.1) and (5.2) the output state is 

(5.3) 
where 

(5.4) ( , 1 = 


Pg (^y )X;^(A'^A + R)^ 


y 

-x' 




Now we specialize the values of u and select the 2n input states 

(5.5) 

(5.6) 

Then the corresponding output states are 


Pj = IV’(2 nej))(i/>(2 1 < j < n, 

Pj = l^(2“^ej))(i/>(2-5ej)|, 1 < j < n. 


(5.7) 
with 

(5.8) 
and 

(5.9) 
with 

(5.10) 


Pg 


-mi 


Pg 


-m; 


lj,mj,-(A'^A + R) 


= A^ 


0 


1 < j < n 




= A' 


0 


1 < j < n. 


A full tomography on (5.7) with j = 1 as outlined in Section (Q yields the hrst row of A and 
the matrix A^A-I-R by using ?7,(2?7,-|-3) measurements. A similar but partial tomography of the 
remaining [n — 1) states in (5.7) and all the states in (5.9) but only for the mean values yields 
the remaining (2?7, — 1) rows of A. This needs an additional set of (2n — l)(2n -|- 1) = 4n^ — 1 
measurements. In all, our approach requires -|- 3n — 1 measurements for getting the QrA + n 
parameters. It will be interesting to know whether one can determine A and B with less 
measurements. If this is not possible our problem will carry an intrinsic tomographic complexity. 
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6. Conclusions 

All the 2n^ + 3n mean and covariance parameters of an n-mode Gaussian state can be recov¬ 
ered from the expectation values of the same number of conjugates of the total number operator 
by Gaussian symmetries. Such symmetries can be realised by hve one mode and four two mode 
gates. The complete tomography of a Gaussian state can be expressed by circuit diagrams and 
measurements akin to those in quantum computation theory. An application of this tomogra¬ 
phy to the output of an u-mode Gaussian channel corresponding to appropriate coherent inputs 
determines all the Qv? + n parameters by 6n^ -|- 3n — 1 measurements. Improvement in this 
channel tomography and hnding the probability distribution of the number operator from its 
explicitly computable probability generating function in a general n-mode Gaussian state seem 
to be interesting problems arising from our investigations. 
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